Uniform semiclassical wave function for coherent 2D electron flow 
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We find a uniform semiclassical (SC) wave function describing coiierent branched flow through 
a two-dimensional electron gas (2DEG), a phenomenon recently discovered by direct imaging of 
the current using scanned probed microscopy [Q, |^. The formation of branches has been explained 
by classical arguments but the SC simulations necessary to account for the coherence are made 
difficult by the proliferation of catastrophes in the phase space. In this paper, expansion in terms 
of "replacement manifolds" is used to find a uniform SC wave function for a cusp singularity. The 
method is then generalized and applied to calculate uniform wave functions for a quantum-map 
model of coherent flow through a 2DEG. Finally, the quantum-map approximation is dropped and 
the method is shown to work for a continuous-time model as well. 



I. INTRODUCTION 

There is no doubt that detailed understanding of the 
electron transport through mesoscopic devices is needed 
to take the full advantage of the possibilities of novel 
electronics these systems offer. On the experimental 
side, great progress was made with the use of scanned 
probe microscopes ^, ^. The theory has kept up: the 
present knowledge has already been summarized in sev- 
eral monographs 1^, ^, 0]. 

Quantum effects have become central as devices have 
become smaller, cooler, and containing fewer impuri- 
ties. Remarkably many quantum properties of the elec- 
tron flow through nanostructures can be explained by SC 
methods. These methods are based on classical mechan- 
ics: the relevant classical manifolds form the "skeleton" 
to which the wave function is attached]^. SC methods 
need to be substituted for classical ones when coherence 
is maintained over distances on the order of the size of 
the device, and when interference effects are playing a 
role. 

In their simplest form, the SC techniques fail when 
nonlinear classical dynamics create complicated struc- 
tures in phase space. In particular, the SC approxima- 
tion breaks down whenever there are multiple contribu- 
tions to the wave function within the volume of a sin- 
gle Planck cell. These so-called catastrophes have been 
classified!^, and various methods have been devised 
to correct the SC wave functions in cases when there ex- 
ist only several coalescing contributions In 
the setting of mesoscopic devices, improved SC methods 
have been applied e.g. to the scattering through ballistic 
microstructures [p^ or to the magnetotransport through 
a resonant tunneling diode |p^ . 

In a recent paper p6[ , we successfully explored a new 
approach which worked even in situations with an infinite 
number of coalescing contributions, occurring e.g. in the 
case of the homoclinic tangle near an unstable periodic 
orbit . This new method is based on the idea of replac- 
ing a complicated classical manifold by a series of new 
simpler manifolds. When standard semiclassical meth- 
ods are applied to these "replacement manifolds" (RMs) , 



accurate uniform wave functions are obtained in situa- 
tions where direct semiclassical evaluation of the original 
manifold fails miserably. 

Originally, this method was used in special, although 
common cases with an infinite number of oscillations with 
the same phase-space area. Here we demonstrate that 
this special property is not necessary, and that similar 
approach may be used more generally, even in cases with 
localized perturbations. In Section II, we briefly review 
the RM method from Ref. and generalize it. The 
method is used to uniformize a cusp singularity in Sec- 
In Section IV , we apply the generalized method 
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tion 

to find a uniform wave function in a quantum-map model 
of a 2D electron flow through a sample with impuri- 
ties, where multiple cusp catastrophes are present. The 
quantum-map approximation is relaxed in Section ^ and 
it is shown how the replacement manifolds are formed in 
a continuous-time model. In Section VI, we discuss the 



merits of the RM method and relate it to other SC tech- 
niques. Because most of this paper is concerned with 
what happens to the twisted manifold under the shear 
of phase space, for completeness the Appendix addresses 
the other major phase space motion: rotation. 



II. REPLACEMENT-MANIFOLD METHOD 
AND ITS GENERALIZATION 

The original method, discussed in detail in Ref. [ p^ , 
works for wave functions of the form 



with 



S{q) = Soiq) + sin /(g). 



(1) 



(2) 



that can be associated with classical manifolds in which 
the momentum depends on the position as 

p{q) ^ ^ + he f{q) cos f{q). (3) 
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Here 5*0 and S are the unperturbed and full action, re- 
spectively, e is a parameter controlling the strength of the 
perturbation, A{q) gives the local weight of the manifold, 
and f(q) is a smooth function defining the shape of the 
perturbation. 

We can expand the wave function as 



n— — oo 



:S„{q) 



(4) 



and interpret each term of the sum as a contribu- 
tion from a classical "replacement" manifold Pn{q) = 
dSn/dq with a weight A„ — A(q)Jn{e) and an action 
Sn{q) = So{q) + nhf{q). The advantage of the RM ex- 
pansion is appreciated after moving to the momentum 
representation with caustics where semiclassical form 
Y,Aj'"'{p)exp[iSj'^'{p)/h] fails while the sum over RMs 
gives an accurate result. 

A slightly different and more general approach than in 
Ref. does not require an oscillatory behavior of the 
action. If 



S{q)^S„{q) + eAS{q), 
we may Taylor expand the wave function as 



(5) 



ipiq) = A(q)exp 



-Soiq) 



oo ^ 
n=0 



h 



eAS{q) 



^An{q) exp 



n=0 



(6) 



corresponding to RMs with weights 

A„(g)=A(<z)(je)"/n! 

and actions 

Sn{q) = So{q) - inn\og[AS{q)/h]. 

Defining a new function f{q) by AS{q) = hexp[f{q)], the 
nth RM action becomes 



Sn{q) = Soiq) ~ ihnf{q) 



(7) 



It will "help" the convergence of expansion (g) if 
limq_>-|-oo fiq) — — oo. This, however, is a natural prop- 
erty of localized perturbations. 

The simplest nontrivial example is obtained by choos- 
ing f{q) = —q^. Besides allowing an analytic solution, 
this choice will yield exactly the manifold needed in our 
model of a 2D electron flow in Section Expanding 
the function p{q) around (7 = 0, 



p{q) = ~2ehqe 



2e%3_q) + 0(q4)^ (g) 



we find that this case falls into the second simplest 
universality class (called cusp) of catastrophe theory 
|, 0, @ (see Fig. 0). 
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FIG. 1: Initial manifold and areas important to the semiclas- 
sical approximation. 



Assuming that the weighting of this manifold is A{q) = 
const. = (27r?i)~^/^, the corresponding SC wave function 
is 



''Psciq) ~ {2TTh) ^1"^ exp (iee 



(9) 



Since for all positions g, there exists only a single contri- 
bution to -^sc (?) ; the SC position wave function is accu- 
rate, ~ ''Psciq)! and the momentum wave function 
is given by the Fourier transform 

i>{p) = {2TTh)-'^^ J dqijsc{q)e-'P'^/^ (10) 

Evaluating this integral by the stationary-phase (SP) 
approximation yields the SC momentum wave function 
ipscip)- The SC momentum wave function has two con- 
tributions from two SP points (Fig. |l|) . The horizontally 
filled- in area gives the phase between two contributions; if 
it becomes smaller than h, the SP approximation breaks 
down. Therefore tpscip) will be singular for all classically 
allowed momenta when e < h. 

Note that the RM momentum p„ (q) — 2inhq is purely 
imaginary for all q and that the corresponding manifold 
has no caustics. The uniform momentum wave function 
is found as 



RM 



/oo 
n=0 



i:S„(g)/ft 



(11) 



S{p) 



2TTh ^ n 

n=l 



-n ' exp 



Sip) + A„e 



where 



Ar. 



2TTh n\ 



-1/2 



Sn{p) = - / dpqn{p) = 



(12) 

(13) 
(14) 



We can in general evaluate all RMs for n > 1 by the SP 
method, although in this case the answer turns out to be 
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equal to the exact Fourier transform because the action 
Sn is quadratic. 



III. UNIFORMIZATION OF A CUSP 
SINGULARITY 

The formation of manifolds with a double-loop struc- 
ture like the one in Fig. |l| is a generic feature of nonlinear 
Hamiltonian systems. This pattern forms, for instance, 
whenever an ensemble of trajectories encounters a dip 
or a bump in the potential surface. Assuming that the 
particles have energy greater the maximum of the poten- 
tial, the dip or bump act as a convex or concave lens, 
respectively. After it is created, the double loop does 
not remain stationary: depending on the Hamiltonian, 
the structure will generally start to shear and rotate in 
phase space (see Fig. ^). In most of this paper we are 
concerned with the shear only, but for completeness, in 
the Appendix we present analytic formulae for the RM 
expansion of an original manifold that is arbitrarily ro- 
tated with respect to the q and p axes. 
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FIG. 2: An example of a manifold with the double loop 
structure that has been sheared (a-c) and rotated (c-d) in 
phase space. While shear and rotation are generic phase 
space motions, here they were implemented by H 
and H = p'^ /2 + /2, respectively. 
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For now imagine that after the manifold (|^) with two 
loops has been formed, the system evolves freely (with 
Hamiltonian H = /2m). Hamilton's equations of mo- 
tion are 



P ■ 



P_ 

m ' 
0, 



(15) 



resulting in a shear of phase-space. The SC position wave 
function . which was accurate at time t = 0, will break 



tc 



m 
2di 



(16) 



when a cusp singularity |l^ develops (see Fig. 

This will be remedied if we apply any of the SC evo- 
lution methods (i.e. integration using the SP approxi- 
mation) to the first few RMs instead of directly to the 
original manifold. 



i'RM{q,t) = J dq'Kf{q,q';t)^RM{q',0) (17) 

oo „ 

= / dq'Kf{q,q';t)A,e^''-^^y^ 



where the free-space propagator 

1/2 



exp 



2ht 



{q" q'f 



and at t = 0, using expression 



-1/2 



{^ey 



ri=0 



(18) 



(19) 



In our case, since the RM terms are Gaussian wave pack- 
ets, their SC evolution (i.e. SP integration) can be per- 
formed analytically and is exact. 



'4>RM{q,t) = (27rft)"'/' 

y ^ (1 + 2inht /to) exp 



(20) 



n=0 



-nq 



1 + 2inht/m 



For comparison, the exact quantum evolution was per- 
formed by switching to the momentum representation, 
using the fast Fourier transform (FFT) and trivially 
evolving the wave function there. To find the primitive 
SC evolution, we used a method described by Berry et 
al.jl^. All three methods are compared in Fig. ||, show- 
ing the classical manifold and corresponding exact, SC, 
and RM wave function at a time instant before, at, and 
after the cusp. 

In the next Section, we show that the RM method 
can treat situations in which more cusps are continuously 
formed. However, the advantage of the RMs over the Van 
Vleck propagation or other standard SC methods can be 
appreciated already when the rough region of the poten- 
tial is localized in time and only one or a few cusps are 
created. As can be seen from Fig. ^, even if the potential 
is simply flat after certain time, the region of q in which 
the SC approximation breaks down expands. Unlike the 
simple SC approximation which deteriorates with time, 
the accuracy of the RM method is preserved after leav- 
ing the rough area of the potential: once the Gaussian 
wavepackets corresponding to the RMs are formed, their 
number remains constant and their propagation is exact 
in any potential with up to quadratic terms (see Fig. |^). 
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FIG. 3: Evolution of the manifold and the comparison of the 
exact (points), RM (solid line) and SC (dashed line) wave 
functions at a time instant before (t — 0.25tcusp, top), at 
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middle), and after {t = Stcusp, bottom) the cusp. 



In these plots, 



1 and the first five RMs were used. 



IV. QUANTUM-MAP MODEL OF A 2D 
ELECTRON FLOW THROUGH A SAMPLE 
WITH IMPURITIES 

We are now prepared to address the problem of the 2D 
electron flow in a semiconductor nanostructure with im- 
purities. The electron transport in such a system is nei- 
ther strictly ballistic nor strictly diffusive. Instead, the 
experiment has revealed that reality lies somewhere in be- 



tween and the phenomenon has been termed "branched 
flow §." Figures I and I show respectively the exact 
electron density (obtained by the exact quantum evolu- 
tion using the FFT) and the representative classical tra- 
jectories in the model described below. The name of the 
phenomenon comes from the shape of the regions with 
enhanced electron density in Fig. ^ or the corresponding 
clusters of classical electron trajectories in Fig. ^j. It turns 
out, however, that these do not correspond to the valleys 
in the potential Although the branches can already 
be seen in the classical simulations. Fig. || also shows that 
scattering by impurities leads to abundant cusp singular- 
ities in phase space, and therefore we expect deviations 
in both classical and primitive SC approximations from 
the exact quantum dynamics. 




t 



FIG. 4: Electron density |i/i(g, t)p in the model of a 2D elec- 
tron flow (obtained by the exact quantum evolution using the 
FFT, RMs were not used). In this plot, e = 2.22 and there 
were 256 impurities. 

Here we analyze a simple model which can neverthe- 
less exhibit all these properties. Namely, we discuss a 
2D system with fast electrons, incident along the x-axis 
and scattered by small isolated Gaussian impurities ran- 
domly distributed in the xy plane. Following Topinka et 
al. 0, 1^ who observed branched flow in a similar sys- 
tem, we consider the electron kinetic energy to be much 
larger than the amplitude of impurities. In fact, we as- 
sume the kinetic energy to be high enough to justify an 
impulse approximation: the electron propagates freely 
between effectively instantaneous kicks from impurities 
that affect its momentum but not position. Moreover, 
while the transverse momentum changes by a small im- 
pulse from an impurity, the longitudinal momentum re- 
mains effectively constant, allowing the transformation of 
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In this section we consider that the change of the trans- 
verse momentum due to the impurity is instantaneous, 
yielding a further simphfication, represented by a peri- 
odically "kicked" Hamiltonian, 




^ t 

FIG. 5; Representative classical electron trajectories, corre- 
sponding to the electron density in Fig. ^. 



the original 2D problem into a ID problem with a time- 
dependent Hamiltonian. To be precise, we start with a 
2D Hamiltonian 



pI 



pI 



E 

i=i 



Vo exp 



2m 



(21) 



(y_y0-))2 



where n, a, and x^-'-' , y^^^ are respectively the number, 
radius, and coordinates of the centers of the impurities. 
We assume that 



Px/m w const. 



(22) 



where v is the initial velocity of the electron. This is 
justified when 



m » v^- 



(23) 



For simplicity of calculations, we distribute the impuri- 
ties randomly in the y direction, but regularly along the 
X axis, at intervals vt. Each electron will be affected by 
a single impurity at a time if 

VT > a. (24) 

To simplify notation, we take a, r, and m to be re- 
spectively the units of length, time, and mass. Defining 
q = y/a, our effective, ID time-dependent Hamiltonian 
becomes 



H{q,p,t) = ?^ 



3 = i 



Vocxp 



(25) 



H{q,P,t) ~ y + VTT— 



\ ^ 



e-(?-?"')^<5(i-j-). 



(A generalized analysis without this approximation is 
presented in the next section.) In the impulse approxi- 
mation, classical position q of an electron does not change 
during an interaction with jth impurity [20[| . Classical dy- 
namics may therefore be expressed in terms of a map, 



qj+i = q, 

Pj+i ^ Pj + ^Pilj+i^l''^^^'') 



(26) 



where subscripts denote time in units r and the change 
of momentum is 



Ap{q,q 



dt 



dH 
dq 



(27) 



implying that a single impurity transforms a momentum 
state exactly into the two-loop manifold (||) from Sec- 
tion II. We can read off the loop area from ( p7| ) to be 
V^\Vo\/v. 

In quantum mechanics, another important parameter 
enters: h. Accuracy of the SC approximation will depend 
on how 



131 

vh 



(28) 



compares to 1. In the impulse approximation, the exact 
quantum dynamics is described by a quantum map 



(29) 



where the subscript again denotes time in units t and U 
is the one-step evolution operator 



U = Texp[-- J Hdt 



expl-l|^Fdtj exp(^-l| 



exp 



exp 



I p 
'h'2 



(30) 



The easiest way to evolve a quantum state numerically is 
to use the FFT to switch back and forth between position 
and momentum representations and apply the impulsive 
part of U in g-representation and the kinetic part of U 
in p-representation. 

We now demonstrate that not only do the replacement 
manifolds lack singularities (present in the classical and 
SC analysis), but that they can also correctly reproduce 



all the details of the exact quantum solution. When the 
next impurity is encountered, each wave packet devel- 
ops a loop in its phase-space representation which would 
soon lead to a new cusp singularity. We therefore replace 
it with a series of simpler manifolds, as in Section 
avoiding this problem. 

In our model we exploit the fact that the RM terms are 
Gaussian wave packets, allowing their analytic evaluation 
with only a slight generalization of the calculations in 
Section III, Each term in the RM sum at time j has a 
Gaussian form, 



ce 



aq~bq 



Re& > 0. 



(31) 



After the kinetic propagation, just before next impurity 
is encountered. 



i^j + i{q) = J Kf{q,q';l)^/j,{q')dq' (32) 

_ c fia^/2 + aq-bq^\ 

" (1 + 2i6)i/2 '^''P 1^ n^2i6 )■ 

After receiving an impulse from the (j -I- l)-st impurity. 



V'i+ilg) = Vij+i(g)exp [iee 



(33) 



n=0 



■ exp 



(g- 9(^+1) )2 



Each term in this su m g ives rise to a new Gaussian wave 
packet of the form (glj); which is propagated further in 
the same manner. 

Figure ^ shows a comparison of the exact and RM evo- 
lution for e = 1.11 up to a time when eight impurities are 
encountered. Four RMs are used to replace each incident 
wave packet at each impurity. Although the classical 
manifold (also shown in the figure) has developed many 
structures smaller than ?i, the agreement remains excel- 
lent. 



V. CONTINUOUS VERSION OF THE MODEL 

In certain situations, we may be interested in a de- 
tailed evolution of the electron wave function during the 
collision with the impurity, rather than just in the ap- 
pearance of the wave function after the collision. Be- 
low, we present an analytical solution of this problem in 
case that the electrons move slowly enough that the colli- 
sion cannot be considered instantaneous, but fast enough 
that the transverse displacement of the electrons does 
not change significantly during the collision. (For even 
slower electrons, the coupling between the longitudinal 
and transverse motion during the collision would prevent 
us from obtaining closed analytic expressions presented 
below. However, we could still find the replacement man- 
ifolds numerically.) 

To simplify the notation, we consider only a single im- 
purity located at position q = and time t = 0, so 






V 



FIG. 6: Evolution of the manifold (left) and comparison of 
the corresponding exact (points) and RM (solid line) wave 
functions (right). In this plot, e = 1.11 and the first four 
RMs were used. Two new impurities were encountered at 
each time interval between the consecutive rows. 



that the effective ID time-dependent Hamiltonian ( |25| ) 
becomes 

H{q, p, t) = ^ (_^2 _ ^2^2^ (34) 

Assuming that q changes little during the collision we 
find that the momentum change is 

Ap = p{t)-p{-^)= f dt'm^- f dt'^ 

J-00 J-00 m 

= 2qFoe"'' / dt' exp {~v^t'^) 

J — 00 

= x/^Vbw"^ge"«' [1 + erf (vt)] (35) 

At time t = —00, we start with a momentum eigenstate 
with momentum p ~ 0, 



ij^ip.t^ -00) = (5(p), 
or ^{q,t = —00) = {2Trh) 



(36) 



represented by a horizontal line in phase space. As the 
electron wave passes through the impurity, a double loop 
develops in the manifold (curve) representing the wave 
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function. The position representation of the SC wave 
function at time t is (see, e.g. [^) 



dqi 



dqj 



1/2 



exp 



n (Si + S2) 



where qi is the position at time t' = —00 that evolves 
to position qf at time t' ~ t. In our approximation qf 
« qi, the Van Vleck determinant \dqi/dqf \ = 1, which is 
the reason that the SC position wave function remains 
accurate throughout the collision. Si is the action along 
the trajectory of a reference point q = x on the manifold, 



^1 



dt'L[x{t'),x{t'),t'] 



S2 is the reduced action along the evolved manifold at 
time t, 



dq'tPt {q't) 



{Pt{q) is the momentum dependence on position at time 
t). For convenience, we choose x(— 00) = —00, giv- 
ing x{t) = and x{t) — const. = —00. Since V{x — 
—00, t) = 0, also L{x, x,t) = and Si = 0. Finally, since 
Pt=-oc{q) = 0, 



^2 



dq' Aptiq') = V^VoV-^ [1 + erf {vt)] 
dq' q' e"«" = Vt^Vov-^^ [1 + erf (vt)] e^'^fs?) 



The semiclassical position wave function at time t is 

(38) 



X exp \ ^VnVov [1 + erf {vt)] e 



Remembering that y/nVov^^h^^ = e and that 
erf (±00) — ±00, we can easily check that this general 
expression gives the correct limiting forms ( ^6| ) and (|^) 
at times t — —00 and t — 00, respectively. The primi- 
tive SC momentum wave function (obtained by the SPA 
of the Fourier transform of (|3^)) fails for the same reasons 
as in Section II. If we expand rpscil^ t) in terms of RMs, 
and apply the SPA directly to the RMs, we find an ac- 
curate answer. The only difference from expression ( |ll| ) 
is an extra factor {| [1 + erf (vt)] }" for RM coefficients 
An or An dTsI), e.g. 



27rft n\ 



/2 J - [1 + erf {vt)] 



Since the expression in the large parentheses goes 
smoothly from at t = — cx3 to 1 at i = 00, we see that 
the replacement manifolds emerge even before the cen- 
ter of the impurity is encountered. However, the weight 
of the manifolds with larger n becomes appreciable only 
after the impurity is passed. 



VI. DISCUSSION AND CONCLUSIONS 



We have shown that the RM method is not limited to 
infinitely repeating phase-space structures if we allow the 
replacement manifolds to have complex momenta. Prop- 
agation of replacement manifolds gives uniform semiclas- 
sical wave functions long after the primitive semiclassical 
approximation breaks down. 

Putting aside the accuracy, the RM approach may 
seem intimidating from a numerical point of view because 
as described, the algorithm has exponential complexity. 
But let us remember that the same-exponential prolifer- 
ation of contributions-is true of the primitive semiclas- 
sical solution-which, however, would give a completely 
wrong result in our case! Moreover, there appear to be 
at least two possible ways to speed up the RM calcula- 
tions: For e < 1 , we could prune the contributions to 
keep only terms up to a certain "total" power of e (which 
is different from keeping all terms up to a given power at 
each impurity). Or we could consolidate the number of 
wavepackets after certain time by projecting on a suit- 
able basis (because the exponentially growing number of 
RM terms is obviously over-complete) and starting the 
RM propagation afresh. 

The question of computational complexity would not 
even arise if we were interested in a system where the 
electron wave hits only one or a few impurities and after 
that propagates in a relatively smooth potential. The 
small number of Gaussian wavepackets spawned at the 
last impurity would suffice for all subsequent times and 
the accuracy of the approximation would be preserved. 
As discussed in Section III , this should be contrasted with 



the standard SC approximation which deteriorates even 
when a manifold with a single cusp propagates in a flat 
potential (see Fig. ||). 

Besides providing a uniform wave function the RM 
method gives an intuitive explanation of how quantum 
mechanics smooths out the classical detail. Moreover, 
in the present case of RMs with a complex momen- 
tum, the method appears to provide a link between the 
semiclassical perturbation approximatio ns |21| and vari- 
ous Gaussian wavepacket techniques!^, |23l , P 4| | , because 
replacement manifolds in the expansion (|L9^are noth- 
ing but Gaussian wavepackets. One advantage of the 
RM method over other Gaussian wavepacket methods 
lies in that it gives an analytic expression for the coef- 
ficients of the wavepackets. Other Gaussian wavepacket 
methods (such as the frozen Gaussians |Q , the Herman- 
Kluk propagator |^ and the full multiple spawning [|4j ) 
rely on variational or ad hoc methods to obtain optimal 
wavepacket coefficients numerically. 

Finally, although the RM method has not yet been 
fully generalized, the large variety of problems (in this 
paper and in Ref. iQ) it can solve suggests that the 
method (or at least, the idea) is more general. 
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APPENDIX A: REPLACEMENT-MANIFOLD 
EXPANSION IN AN ARBITRARILY ROTATED 
COORDINATE SYSTEM 

We show here that the RM expansion for the manifold 
studied in this paper can be found analytically in an ar- 
bitrarily rotated coordinate system. In other words, the 
method can be readily applied not just in the q or p rep- 
resentations, but in any mixed representation given by a 
canonical transformation (see Fig. 0) 



Q ^ q cos (f) + p sin (f) 
P = —q sin (j) + p cos ( 



(Al) 



n=2 




FIG. 7: Original and rotated coordinate systems. 



In the original coordinate system {q, p), replacement 
manifolds are straight lines 



Pn{q) — 2inhq = q tan an 



(A2) 



(where is complex). In the rotated coordinates (Q, 
P), the replacement manifolds are given by the relation- 
ship 



PniQ) = (5tan(a„ 
The reduced action is 



Q 



2inh — tan(/) 
1 -|- 2m/itan( 



(A3) 



SniQ) - J dQPniQ) = ^QHanian 



0). (A4) 



The weight of the nth RM in the Q representation is the 
weight in the q representation multiplied by the ratios of 
the projections on the q and Q axes, respectively. Includ- 
ing the Maslov index /z (0 or 1), we find the correct nth 
RM contribution 



i>RM,n{Q) 



1 



'2'Kh n 



cos 



cos(a„ - 0) 



1/2 



X exp 



2h 



' tan(a„ — 0) — i/i„7r/2 



(A5) 



After simplification, the full RM expansion becomes 



iter 



(A6) 



E 



/2Trh ^ nl 

n—O 



I cos -I- 2in/isin</)| 



X exp 



,n + i tan (f)/(2h) 
1 + 2m?itan ch 



i^„7r/2 



It is easy to check that this general result correctly re- 
duces to expressions (0) or (|l|) when (/> = 7r/2 or = 0, 
respectively. 
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